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ABSTRACT 

Previous surveys of outflows in low-redshift active galactic nuclei (AGN) have relied 
on the analysis of sources selected primarily for their optical/X-ray brightness, and are 
therefore biased. Towards determining the outflow properties of local AGN, we detect 
warm absorption signatures of O VII and O VIII absorption edges in the available 
Suzaku/XMM-Newton CCD spectra of an unbiased sample of 44 Seyfert 1-1.5 sources 
selected in the very hard X-rays (14-195 keV) with the Swift Burst Alert Telescope. 
From our analysis, we find that O VII and O VIII absorption edges are present in 41% 
of the sample. This fraction is dependent on luminosity, with outflow detections in 60% 
of low luminosity and 30% of high luminosity sources. However, grating spectroscopy 
of the highest luminosity sources reveals that ~ 80% of these sources have ionized 
absorbers, but that the ionization states are higher/lower than produces the O VII and 
O VIII edges. This suggests that ionized absorption may be present in all local Seyfert 
Is. 

Subject headings: galaxies: active — galaxies: Seyfert — X-rays: galaxies 



1. Introduction 

Mass outflows from active galactic nucleus sources (AGNs) are believed to affect their sur- 
rounding galaxy and the intergalactic environment in many ways. The powerful winds from AGN 
may play an active role in shaping the host galaxy evolution by blowing out the necessary fuel 
for star formation, effectively quenching it. Further, the ejected gas enriches the galaxy's local 
environment as well as the intergalactic medium with metals produced from supernovae. While 
it is believed that outflows are ubiquitous in AGN, signatures of t hese winds are only seen in 



50-60% of Seyfert (Sy) 1—1 .5 sources in the X-ray and UV bands (jRevnolds! 1 19971 : iGeorge et al 



1998 



Crenshaw et al 



19991 ) . This fraction is of fundamental importance, since it represents the 



solid angle occupied by the outflow. 

The previous outflow studies are not complete - selecting bright sources in the optical/soft X- 
ray bands. Therefore, they may be biased towards selecting warm absorbers. To determine the true 
fraction of AGN with outflows/warm absorbers, it is necessary to search for absorption signatures 
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in an unbiased sample. In this letter, we detail our study identifying outflows in the sources selected 
in the very hard X-ray band with Swift's Burst Alert Telescope (BAT). Due to the selection in the 
14-195 keV band, the Swift sample is unbiased towards all but the heaviest line-of-sight absorption 
(Nh > 10 24 atoms cm -2 ) for sources with a BAT flux of > 2 x 10 -11 ergss -1 cm -2 . 



2. The Swift Sample 



The sample of 153 AGNs in the Swift BAT nine-month catalog (jTueller et al.l 120081 ) are 
the most well-studied a nd brightest hard X-ray sources. De tails of their 0.3-10 keV X-ray spec- 
tra (IWinter et al.l 120091 ). optical spectra (jWinter et al.ll2010l ). and infrared spectra from Spitzer 
(jWeaver et al.ll2010l ) are well documented in the literature. Additionally, many of these sources are 
well-known AGNs with high signal-to-noise archived data available from current X-ray observatories 
(i.e., Suzaku, XMM-Newton, Chandra). 

We selected all of the optically classified Seyfert 1—1.5 sources in the nine-month BAT catalog 
with a high Galactic latitude (|6| > 15°). This sample includes 51 sources, of which 44 have archived 
spectra available. The typically low intrinsic column densities (Nh < 10 21 cm -2 ) of Seyfert Is make 
them the best sources for searching for outflow signatures. This is due to the fact that X-ray outflow 
signatures are found primarily in the soft band, where more absorbed sources have lower observed 
luminosities. 



3. X-ray Observations 

In the X-ray band, outflows are detected primarily in the soft 0.3-2 keV band through emission 
and/or absorption features from oxygen (e.g., 0.65 keV O VIII Lya, 0.77 keV O VIII Ly/3), magne- 
sium (e.g., 1.35 keV Mg XI), neon (e.g., 1.02 keV Ne X Lya), and iron (e.g., 0.87 keV Fe XVIII, 0.92 
keV Fe XIX, 0.96 keV Fe XX). However, detection of these lines typically requires high signal-to- 
noise observations with the grating spectrometers on XMM-Newton or Chandra. To obtain high 
signal-to-noise grating spectra, long exposures are required of the order of > 100 ks for the Swift 
sources. Therefore, it is more feasible to detect X-ray outflows in a larger sample of sources through 
CCD data. In the CCD data, absorption edges from O VII (0.73 keV) and O VIII (0.87 keV) indicate 
the presence of outflows. 

High signal-to-noise archived X-ray CCD data from Suzaku and XMM-Newton are available 
for many of our sources. While the O VII and O VIII edges are d etectable in data o f moderate 



quality (i.e., ~ 10 ks with XMM-Newton for the Swift sources as in lWinter et al.ll2008l ). we chose 
to use predominantly higher exposure time observations of > 40 ks (though, exposure times range 
from 2-275 ks with 34/44 having exposure times > 40 ks). The Suzaku X-ray observatory was our 
first choice for X-ray observations since simultaneous data is available from the 0.3-12 keV (with 
the XIS detectors) and ~ 15-50 keV band (with the HXD PIN detector). This allows us to better 
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constrain the underlying continuum and provides overlapping energ y coverage between the Suzaku 



spectra and the 22-month averaged 14-195 keV Swift BAT spectra (jTueller et al.ll2010l ). Properly 
constraining the continuum emission is essential in order to accurately determine limits on the 
optical depth of O VII and O VIII absorption edges. Suzaku data were available for 27/51 sources. 
Where multiple Suzaku observations exist, we analyzed the longest observation. 

For sources without a Suzaku follow-up, we analyzed archived XMM- Newton EPIC observa- 
tions. An additional 17 sources have high signal-to- noise XMM-Newton observations without a 
Suzaku follow-up. In total, we analyzed the spectra of 44/51 (80%) of the Seyfert 1-1.5 sources 
in the Swift BAT nine-month catalog. The sources which do not have available high quality ob- 
servations represent lower luminosity sources which were identified as AGN for the first time by 
Swift. 

Full details of the observations and data are presented in Winter et al. (in prep). To summa- 
rize the processing of the observations, for the Suzaku XIS (XISO, XIS1, XIS3, and XIS2, where 
available - in the 0.3-12 keV band) and XMM-Newton EPIC pn (in the 0.3-10 keV band) data 
we extracted the spectra from circular regions centered on the source of ~ 60-200". Background 
regions were extracted in regions near the source but free of additional sources. The Suzaku spectra 
were extracted using the ftool xselect, while the pn spectra were extracted with SAS using the 
standard filtering for the pn as detailed in the SAS ABC Guide. Response and ancillary response 
matrices were generated using the Suzaku ftools xisrmfgen and xisarfgen or SAS commands 
rmfgen and arfgen. The spectra were then grouped and binned to 20 counts per bin using grppha. 



The Suzaku HXD PIN (~ 15-50 keV; iTakahashi et al.l 120071 ) spectra were extracted using 



xselect from the processed cleaned event files. The tuned PIN background file specific to each 
observation was used to extract a background spectrum. The spectra were corrected for dead time 
and the CXB contribution was added to the background spectrum as specified in the Suzaku ABC 
Guide. The corresponding response file from the Suzaku CALDB was used for each observation. 
The PIN spectra and responses were grouped using grppha and binned to a signal-to-noise of 3-4<r. 
Finally, we also util ized the Swift BAT 22-month averaged spectra and diagonal response file from 



Tueller et al.l (J2010J). These spectra consist of 8 energy channels in the BAT band of 14-195 keV. 



4. Spectral Analysis 

We simultaneously fit the Suzaku XIS and PIN or XMM-Newton pn spectra with the Swift 



BAT spectra to obtain joint fits in the 0.3-195 keV band. We used the XSPEC (jArnaudl 1 19961 ) 
spectral fitting package vl2. For the sources with Suzaku spectra, we added a constant model to 
allow for differences in the flux scaling between the spectra, normalizing to the XIS1 spectrum. 
The PIN spectrum was fixed to the value of 1.16, as detailed in the Suzaku ABC Guide, while the 
constant factor was allowed to vary for the BAT spectrum and the additional XIS spectra. For 
sources with XMM-Newton spectra, we allowed the constant factor to vary relative to the pn level. 
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For all of the sources, we fit the broad-band X-ray spectra with a base cutoff power law model. 
Seyfert direct emission is assumed to be produced from inverse Compton scattering in the corona 
surrounding an accretion disk - which is well represented by a power law with a cutoff at high 
energy. We included n eutral absorption from the Galaxy with a tbabs model, fixing the Galactic 



column density to the iDickev Lockmanl ([1990) value. Where the addition of intrinsic neutral 



absorption improves the fit (Ax 2 > 6.64 on adding the model, corresponding to P = 0.01 for the 
addition of one degree of freedom), we added intrinsic emission with a ztbabs model. Where a 
soft excess was evident, we included a blackbody model to account for the soft excess. To better 
constrain the continuum emission, we added a Gaussian model (zgauss) to account for neutral Fe 
Ka emission and a pexrav model to account for reflected emission from the accretion disk. Full 
details of the modeling, including the best-fit parameters (in particular, L 14 _ 1 9 5 k e v and Nh), are 
in Winter et al. (in prep). 

With a baseline model for the continuum, we searched the spectra for signatures of outflows 
through the presence of O VII and O VIII edges. We used two zedge models, one to account for 
each of the two absorption edges. The zedge model has three parameters, the energy of the edge, 
redshift, and optical depth. We initially allowed the edge energy and optical depth to vary around 
the laboratory measured values of the edges (0.73 and 0.87 keV, respectively), fixing the redshift 
to the source redshift. Where the energies were not well-constrained (i.e., the fitted energies were 
not consistent with the absorption edge energies), we fixed the energies to the laboratory value. 
In Table [IJ we record the results of this analysis, including the A% 2 upon adding the absorption 
edges to the baseline model, energies, and optical depths for the fitted edges. We classify the source 
as exhibiting an outflow where an improvement in x 2 corresponds to Ax 2 13.39 (a probability 
of P = 0.01 for four additional degrees of freedom). Additionally, we include NGC 3516, whose 



spect rum was not well-fit by our base model, but clearly has an outflow present (see lMarkowitz et al 



20081 for a complete analysis of the absorption in the Suzaku observation of NGC 3516). 



5. Results 

We classify 41% (18/44) of Swift BAT Sy 1-1.5 sources as exhibiting the outflow signatures 
probed through detection of O VII and O VIII edges. Of the outflow sources, half are Sy 1-1. 2s 
and half are Sy 1.5s. The range of optical depths for O VII is t ~ 0.07 — 2.54 and for O VIII is 
r ~ — 1.66. For non-outflow sources, we find the range of optical depths for O VII is r ~ — 0.29 
and for O VIII is r ~ — 0.26. In this section, we describe the properties of the outflow versus 
non-outflow sources. 

In Figure 1, we plot the optical depth of O VII and O VIII versus the Swift BAT luminosity 
(14-195 keV) for the outflow (circle) and non-outflow (triangle) sources in our sample. Throughout 
this section, we use the Swift BAT luminosity because it is measured at high energies and is thus 
less affected by obscuration. Clearly, the outflow sources have the highest measurements for optical 
depth, on average. Additionally, we find that the measured optical depth in both O VII and O VIII 



- 5 - 



is low for the highest luminosity (£i4-i95keV ^ 5 x 10 43 ergss _1 ) sources overall. We plot the Sy 
1-1.2 sources with open and the Sy 1.5 sources with filled symbols. Another trend which we find 
is that the highest luminosity sources tend to be Sy 1-1.2 sources. Only one source out of eleven 
with Li4_i95k e v ^ 10 ergss -1 is a Sy 1.5. Finally, the four sources with the strongest outflows 
(NGC 3516, NGC 4151, Mrk 6, and NGC 526A) have intermediate luminosities and are all Sy 1.5s. 

Additionally in Figure 1, we plot the relationship between optical depth and accretion rate 
for our sample. As an estimate of accretion rate, we use Li4_i95k e v/LEdd> where LEdd = 1-38 x 
10 38 M/M(^ er g ss -1 . The black hole mass estimates are from stellar bulge K-band photometry in 



Winter et al.l (|2009l) , which were show n to be well-correlated with reverberation mapping based 



mass estimates (see lWinter et al.ll2010l ). There is no correlation between accretion rate and optical 
depth. However, the sources with the strongest outflows have intermediate values of accretion rates 
and luminosities. The highest accretion rate sources are Sy 1—1.2 sources. 

Table [2] includes average, standard deviation, and the results of Kolmogrov-Smirnov (K-S) 
statistical tests for measured parameters of our outflow and non-outflow samples. Comparing the 
average values of both luminosity and accretion rate, we find that the sources with outflows have 
lower luminosities and accretion rates than the sources without outflows. Also, we find that the non- 
outflow sources have more massive black holes. There is also a difference between the outflow and 
non-outflow sources in the measured neutral hydrogen column density. Assuming Nh = 10 19 cm~ 2 
for sources with no significant Nh measured, we find that sources with detected outflows have 
larger neutral column densities. 

Results of the K-S tests show that the probability of mass and accretion rates of the outflow and 
non-outflow sources being drawn from different populations is low. Therefore, the main statistical 
difference between the outflow and non-outflow sources is the luminosity and secondarily Nh- 

In Figure 2, we plot the fraction of outflows detected in binned Nh, BAT luminosity, black hole 
mass, and accretion rate. We find that there are few outflows (< 30%) detected in sources with 
low column densities (Nh < 10 20 cm -2 ). Conversely, outflows are detected in ~ 63% of sources 
with Nh > 10 20 cm -2 . Outflow detection rates are low at the highest luminosities (~ 20% for 
Li4_i95keV ^ 10 44 ergss" 1 ), masses (~ 27% for M/M Q > 3 x 10 8 ), and accretion rates (~ 27% for 
Li4_i95kev/LEdd ^ 0.003). Therefore, even though the total distribution of black hole mass and 
accretion rate between outflow and non-outflow sources is not statistically different, the binned 
values show fewer outflows at the highest mass/accretion rate. At the lowest masses/accretion 
rates/luminosities, we find outflows in ~ 50-60% of the sources - approximately twice the detection 
rate for the highest mass/accretion rate/luminosity sources. 



6. Discussion 



We have categorized sources in our study as outflow or non-outflow sources using the signif- 
icance of O VII and O VIII absorption edges to a base continuum model (an absorbed power-law 
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with soft excess and reflection). This ensures that we are selecting the sources with the highest 
optical depth in VII and O VIII as outflow sources (see Figure 1). Our sample includes 44 Sy 
1—1.5 sources detected in the very hard X-rays with the Swift BAT. This is the first study of the 
occurrence of outflows in an unbiased sample of AGN. 

We find an outflow detection rate of 41%. In a similar study detecting outflows thro ugh O VII 
and O VIII edges in ASCA spectra of 24 type 1 AGN, a detection rate of 50% is found ([Reynolds 
19971 ) . Of these 24 AGN, 18 are also found in our survey. We find the same categorization as an 
outflow/non-outflow source for 14/18 sources. The remaining four (IC 4329A, NGC 5548, 3C 382, 
and 3C 390.3) are categorized as having outflows in lRevnoldsl (jl997l ) and non-outflow sources in our 
study. The upper limits on the optical depths for IC 4329A and NGC 5548 are within the range 
of values for sources with detected outflows, though the A% 2 values are low. Further, Chandra 



X-ray grating observations of these sources show ionized absorption is present ([Steenbrugge et al 



2005bl lal). The residuals to our base model show evidence of potential warm absorber signatures in 
the spectra of IC 4329A and NGC 5548, but at lower energies (< 0.7 keV) than probed with the 
O VII and O VIII edges. These features are likely associated with absorption in N VII and O VII 
from ~ 0.50 - 0.55 keV. 

Grating spectroscopy with XMM- Newton of 3C 382 re vealed a highly ionized (log £ = 2.69 erg cm s~ 
absorber with a column of 10 22 cm -2 ( Torresi et al.ll2010i ). Finally, ionized absorption is detected 
at higher energies than probed in our study (6.6 keV) in XMM-Newton o bservations of 3C 390.3 , 
indicating an even higher ionization parameter of log£ = 3.43 erg cms -1 ( Sambruna et al.ll2009l ). 
Therefore, we miss the detection of outflows in these luminous sources because they are more highly 
ionized. It is unclear, however, why O VII and O VIII edges were detected in the ASCA data, but 
not with the Suzaku/XMM-Newton CCD data. 

The main question that our results raise is why there are less outflow detections in the most 
luminous AGN. The most luminous sources in our sample, those with log £14-195 keV ^ 44.5, are, 
from most to least luminous, 4C +74.26, MR 2251-178, 3C 390.3, 3C 382, 1H 0419-577, and Mrk 
926. The four highest luminosity sources are radio-loud and have very broad optical emission line 
profiles, including double-peaked hydrogen Balmer lines. Very highly ionized (log£ > 3) gas is 
detected in each of the se sources ( grating spectroscopy reveals highly ionized gas in 4C +74.26 
and MR 2251-178, see iKaspil 120041 ) . Of the additional luminous sources, XMM-Newton grating 
spectroscopy of 1H 0419-577 reveals ionized gas with a low ionization parameter (lo g£ ~ 1.3) and 



mult iple observations show that the ionized (as well as neutral) absorption is variable (jPounds et al 



2004). Since high resolution observations of Mrk 926 do not exist, it is unclear whether ionized 
absorbers are present in this source. However, it is clear that the fraction of ionized absorbers 
is high even in the highest luminosity sources (at least 5/6), despite the fact that O VII and 
O VIII absorption edges are only detected in 2/6 of these sources. Similarly, 3/4 sources with 
log Li4_i95keV < 42.5 have absorbers detected (two are outflow s ources in this paper , NGC 7213 
has detected OVII and OVIII emission lines in HETG spectra ([Starling et al.ll2005l ). and UGC 
6728 has no high resolution spectra available). 
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The next step in understanding the outflow properties of local Seyfert Is is to identify the more 
detailed properties of the ionized gas (i.e., ionization parameter, column density, velocity). Our 
study of outflow detections from O VII and O VIII absorption edges probes only a narrow range of 
ionization states. As we found for the highest luminosity sources, we are missing the detection of 
more or less ionized gas and therefore underestimating the fraction of sources with outflows. Based 
on a literature search of the highest luminosity sources, we find that the fraction of outflows is 
> 80%, while the fraction of sources with O VII and O VIII absorption edges is from 30-60% (with 
low detection rates at high luminosity and high detection rates at low luminosity). This suggests 
that outflows may be present in all local AGNs. As future work, we will extend our outflow study 
by analyzing archived grating spectroscopy to determine the ionization state, column density, and 
velocity of the warm absorbers present in the Swift BAT sample. 

LMW thanks Tatiana Taylor for assistance with data reduction. She also acknowledges support 
through NASA grant HST-HF-51263.01-A, through a Hubble Fellowship from the Space Telescope 
Science Institute. Facilities: SwiftQ, XMM(), SuzakuQ 
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Table 1. Warm Absorption through the Detection of VII and VIII Absorption Edges 
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Table 1 — Continued 



Source 


Type 


Ax 2 


OVII 

E (keV) 




OVII 

r 


OVIII 

E (keV) 




OVIII 

r 


EXO 055620-3820.2 


1 


.0 


-4.19 


0. 


73 


0. 


•000+0 


.047 
.000 


0. 


.87 







nnn+o- 012 

uuu -0.000 


2MASX J09043699+5536025 


1 


.0 


12.99 


0. 


.73 


0. 


• 267 -a 


.204 
.205 


0. 


on+0.04 
°°-0.04 


0. 


25Q+0-141 


MCG +04-22-042 


1 


.2 


0.21 


0. 


.74 


0. 


.000+°, 


.013 
.000 


0. 


,87 




0. 


m 9 +0.029 
u±z -0.012 


Mrk 110 


1 


.0 


6.68 


0. 


,73 


0. 


.124+0; 


.036 
.040 


0. 


,87 




0. 


onn +o.oii 

uuu -0.000 


UGC 06728 


1 


.2 


3.8 


0. 


.73 





•023±g; 


.058 
.023 


0. 


,87 




0. 


m 7 +0.037 
ul ' -0.017 


IC 4329A 


1 


.2 


-11.69 


0. 


.73 


0. 


,263+°, 


.052 
.054 


0. 


o fi +0.01 
86-0.01 


0. 


1 07+0. 036 
±0 ' -0.023 


Mrk 279 


1 


.5 


1.88 


0. 


.73 





.014+0 


.016 
.014 


0. 


.87 




0. 


nnn +0.004 

• uuu -o.ooo 


NGC 5548 


1 


.5 


-1.96 


0. 


,73 





,054+0 


,105 
.054 


0. 


.87 




0. 


nnn +0.065 
uuu -0.000 


ESO 511-G030 


1 


.0 


-0.29 


0. 


.73 


0. 


.oooi , 


.006 
.000 


0. 


.87 




0. 


nnn +0.002 
uuu -0.000 


Mrk 841 


1 


.0 


13 


0. 


7^+0.04 
'°-0.03 


0. 


.0671°, 


.045 
.045 


0. 


,87 




0. 


nnn +0.030 
uuu -0.000 


3C 382 


1 


.0 


7.46 


0. 


.73 


0. 


.oooi , 


.003 
.000 


0. 


OO+0.03 
OO-0.03 





n2n +o.oio 

uzu -0.010 


3C 390.3 


1 


.0 


-0.32 


0. 


.73 





.oooi , 


.007 
.000 


0. 


.87 







nnn +0.008 

• uuu -o.ooo 


NGC 7213 


1 


.5 


-59.3 


0. 


.73 


0. 


.oooi , 


.005 
.000 


0. 


.87 




0. 


010+0- 016 

Uiu -0.010 


Mrk 926 


1 


.5 


-1.53 


0. 


.73 





.oooi , 


.025 
.000 


0. 


.87 




0. 


nnn +0.022 

• uuu -o.ooo 



This table includes results of fitting zedge models to the O VII and O VIII absorption edges at 

0.73 keV and 0.87 keV, respectively. Type is the optical Seyfert type for each AGN. The Ax 2 
measurement represents the improvement of the fit upon adding the edges to our baseline model. 

We record both the measured edge energies (keV) and optical depths (r). Where the energies 
were not well-constrained, we fixed these to the lab value. For NGC 3516, NGC 4051, and NGC 

4151, our baseline model was not sufficient to obtain a good fit to their spectra (x 2 /dof > 2). 
Therefore, error bars are not computed for these sources. 



Table 2. Comparison of the Properties of Outflow vs. Non-outflow Sources 



Parameter 


Outflow 


Non-outflow 


K-S D 


K-S P 


log L 1 4_ 195ke v 


43.38 ± 0.79 


43.78 ± 0.63 


0.5256 


0.003 


log Li4_i95kev/L.E(i(i 


-2.92 ± 0.49 


-2.80 ±0.64 


0.2265 


0.588 


log M/M 


7.80 ± 0.52 


8.09 ± 0.63 


0.3419 


0.131 


log N H 


20.57 ± 1.12 


19.42 ± 0.99 


0.4487 


0.018 



This table includes the average and standard deviations on indicated parameters for the outflow 
and non-outflow sources. Units for luminosity and Nh are ergss -1 and cm -2 , respectively. Both 
the K-S D parameter and probability for the distributions being drawn from different populations 

are also given. 
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Fig. 1.— Plotted is the optical depth of O VII (left) and O VIII (right) versus the Swift BAT 
luminosity (top) and the ratio of the Swift BAT luminosity to the Eddington luminosity (bottom). 
Sources where A% 2 significantly improved upon adding the edge models (outflow sources) are shown 
with circles. Non-outflow sources are represented by triangles. Additionally, open symbols represent 
Sy 1-1. 2s, while filled symbols represent Sy 1.5s. The highest luminosity sources, as well as sources 
with the highest accretion rates, tend to be Sy 1-1. 2s. 
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Fig. 2. — Plotted, we show the fraction of outflow sources/total sources in binned intrinsic neutral 
column density (top left), Swift BAT luminosity (top right), black hole mass (bottom left), and 
accretion rate (bottom right). We label the number of outflow sources out of the total number 
of sources in each bin. Outflows tend to be detected in sources with higher Nh- There are few 
outflows detected in the most luminous sources, sources with the most massive black holes, and 
sources with the highest accretion rates. 
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